Random walks in modular scale-free networks with multiple traps 
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Extensive empirical investigation has shown that a plethora of real networks synchronously exhibit 
scale-free and modular structure, and it is thus of great importance to uncover the effects of these 
two striking properties on various dynamical processes occurring on such networks. In this paper, 
we examine two cases of random walks performed on a class of modular scale-free networks with 
multiple traps located at several given nodes. We first derive a formula of the mean first-passage time 
(MFPT) for a general network, which is the mean of the expected time to absorption originating from 
a specific node, averaged over all non-trap starting nodes. Although the computation is complex, 
the expression of the formula is exact; moreover, the computational approach and procedure are 
independent of the number and position of the traps. We then determine analytically the MFPT for 
the two random walks being considered. The obtained analytical results are in complete agreement 
with the numerical ones. Our results show that the number and location of traps play an important 
role in the behavior of the MFPT, since for both cases the MFPT grows as a power-law function 
of the number of nodes, but their exponents are quite different. We demonstrate that the root of 
the difference in the behavior of MFPT is attributed to the modular and scale-free topologies of 
the networks. This work can deepen the understanding of diffusion on networks with modular and 
scale-free architecture and motivate relevant studies for random walks running on complex random 
networks with multiple traps. 

PACS numbers: 05.40.Fb, 89.75. He, 05.60.Cd, 89.75.Da 



I. INTRODUCTION 

In the past decade, with a huge amount of data and 
computational resources available, scientists have pro- 
cessed and analyzed data of a wide variety of real sys- 
tems in different areas, leading to important advances in 
the understanding of complex systems [IH1] . A large vol- 
ume of empirical studies showed that scale- free feature [5[ 
and modular structure [6-8] are two prominent proper- 
ties that seem to be common to real networks, especially 
biological and social networks. The former implies that 
the networks obey a power-law degree distribution as 
P(k) ~ fc -7 with 2 < 7 < 3, while the latter means 
that the networks can be divided into groups (modules) , 
within which nodes are more tightly connected with each 
other than with nodes outside. These two remarkable 
natures constitute our fundamental understanding of the 
structure of complex networks, which are relevant to 
other topological features (i.e., average distance [H, [1(3] 
and clustering coefficient Q), and have led to many pop- 
ular topics of research in network science, including ex- 
plaining the origin of the scale- free phenomenon 0, 
identifying the modules [IiT - [l5l ] and finding their con- 
crete applications [H, E3 • 

It is well known that one of the ultimate goals for re- 
search on complex networks is to make clear how the 
underlying structural characteristics affect the dynami- 
cal processes defined on networks [IGJ]. Among various 
dynamical processes, random walks have held continual 
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interest within the scientific community [19-35] because 
of their relevance to a wide range of different applica- 
tions to many fields (36l . |37| . In particular, as an inte- 
gral subject of random walks, the trapping problem is 
closely related to numerous aspects in a great many dis- 
ciplines [38Tj4l| . Over the past decades, scholars in a 
large interdisciplinary community have made a huge ef- 
fort to address the trappin g pr oblem in diverse networks, 
including regular lattices [42j], regular fractals [4 
small-world networks 1481. a nd scale-free networks 
[55| . among other graphs |56l - l58j . 

Thus far, most previous works on random walks in 
complex networks have focused on the case with a sin- 
gle trap fixed at a given location, while work on the case 
with multiple traps is much less common. In particu- 
lar, research on the multiple-trap problem in complex 
networks with modular organization and scale-free struc- 
ture is still lacking, despite the multiple-trap issue having 
obvious applications to various as pec ts [59| (description 
of particle-cluster aggregation (fjfji. l6l|. for instance) and 
being relevant in diffusion-limited reactions in chemical 
field [62| and modular and scale-free topologies having 
vital influence on dynamical processes taking place on 
networks HHfTZj]- 

In this paper, we study the classic random-walk prob- 
lem for a category of modular scale-free networks 0, ||| 
with several given nodes being occupied by immobile 
traps, which absorb all particles visiting it. The basic 
quantity we are interested in is the mean first-passage 
time (MFPT) [68[ characterizing the trapping process, 
which is defined as the average of expected time for a 
particle starting off from a particular node until first vis- 
iting one of the traps, averaged over all nontrap source 
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nodes. The networks we study are of a deterministic fam- 
ily, which has proven to be an important tool in the field 
of complex networks and has recently attracted much in- 
terest (si-iil. 

We first study the MFPT in a generic network and re- 
duce the problem of computing the MFPT to finding the 
sum of elements of a matrix associated with the trapping 
issue. Then, based on the deterministic recursive con- 
struction of the modular scale-free networks being con- 
sidered, we investigate analytically the key of quantity 
MFPT for two cases of particular arrangements of traps. 
In the first case, traps are placed on peripheral nodes; 
in the other case, traps are fixed on those nodes far- 
thest from the hub. For both cases, we derive exactly 
the dominant scalings for the MFPT and show that they 
produce a power-law function of the network size with 
their exponents smaller than 1 but being different, which 
is confirmed by the numerical results obtained via invert- 
ing related matrices. The obtained results indicate that 
both the number and the location of traps have a signifi- 
cant impact on the behavior of the trapping. We demon- 
strate that the high efficiency of both trapping processes 
and the distinction between the behavior of the MFPT 
for the two random walks are rested with the scale-free 
property and the modular structure of the networks un- 
der consideration. 



II. MODULAR SCALE-FREE NETWORKS 




FIG. 1: (Color online) Sketch of a network M3 for the limiting 
case of m = 5. Note that the diagonal nodes are also linked 
to each other; the edges are not visible. 

tinct sets [H|, HH : the peripheral node set P, the locally 
peripheral node set P z (1 < z < g), the set H contain- 
ing only the hub node of M g , and the local hub set H z 
(1 < z < g). The number of nodes in each of these four 
sets is 



I IP I — ( m - 1)9 (1) 

We first introduce the model for the modular scale-free 
networks, which are built in an iterative way Let 

M g stand for the network model after g {g > 1) iterations |p | _ r m _ iy m 9-( z + 1 ) (2) 

(i.e., number of generations) . Initially (g = 1), the model 
is composed of to (to > 3) nodes linked by m(m — l)/2 

edges forming a complete graph, among which a node |H| = 1, (3) 

(e.g., the central node in Fig. [1]) is called hub (or root) 

node, and the other to — 1 nodes are named peripheral and 

nodes. At the second generation (g = 2), m — 1 repli- ^ (z+i) 

cas of M\ are created with the m — 1 peripheral nodes 1^1 = ( m ~ l) m9 , (4) 

of each copy being connected to the root of the original 

Ml In this way, we obtain M 2 , the hub and peripheral respectively. For M g , all nodes m a set have the same 
nodes of which are the hub of the original M 1 and the de S ree ' Tt has been obtained exactly that the degree for 
(to - l) 2 peripheral nodes in the m - 1 duplicates of M x , a nodc 111 scts H ' H *' P ' and P * 1S ' respectively, 
respectively. Supposing one has M g _i, the next genera- g 

tion network M g can be obtained by adding to — 1 copies Kh{g) = (to — \) 9i = m — [(m — I) 9 — 1] , (5) 

of Mg_i to the primal M g _i, with all peripheral nodes of g . =1 m — 2 

the replicas being linked to the hub of the original M g _i 
unit. The hub of the original M g _\ and the peripheral 

nodes of the to — 1 copies of Af g _i form the hub node and Tr / N / „„. m — 1, 

peripheral nodes of M g , respectively. Repeating indefi- 
nitely the two steps of replication and connection, one 
obtains the modular scale-free networks. Figure Q] illus- 
trates a network M4 for the particular case of m = 5. 

Many interesting quantities and properties of the 
model can be determined explicitly [H, [82|. In M g , an d 
the network size (number of nodes), denoted by N g , is 

N g = to 9 . All these nodes can be classified into four dis- K PtZ (g) = z + m — 2. (8) 



K h ,M = E ( m - l ) m = ^[(™ - !) Z - 1] . ( 6 ) 

^—^ Tfl — 2 

K P (g) =g + m-2, (7) 
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In addition, it is easy to obtain that the average degree 
of all nodes is approximately equal to a constant 2(m — 
l)(3m — 2)/m in the limit of infinite g, showing that the 
networks are sparse. 

The model under consideration is in fact an extension 
of the one proposed in Ref. [HJ and studied in much 
detail in Refs. |85l - l87| . It presents some typical features 
observed in a variety of real- world systems [H, HH . Its de- 
gree distribution follows a power-law scaling P(k) ~ k 1 
with a general exponent 7=1 + lnm/ln(m — 1) be- 
longing to the interval (2,2.585]. Its average clustering 
coefficient tends to a large constant dependent on m; and 
its average distance grows logarithmically with the net- 
work order I88I1. both of which show that the model is 
small world [89(. In addition, the betweenness distri- 
bution P(b) of nodes also obeys a power-law behavior 
P(b) ~ 6~ 2 with the exponent independent of the pa- 
rameter m. In particular, the whole class of networks 
shows a remarkable modular structure. These peculiar 
structural properties make the networks unique within 
the category of complex networks. It is thus interesting 
to address dynamical processes happening on them. The 
main purpose of this work is to study random walks on 
this network family with multiple traps located on some 
special nodes. 



III. FORMULATION OF RANDOM WALKS 
WITH MULTIPLE TRAPS ON A NETWORK 



It is clear that for any node i £ T, we have Tj = 0. The 
set of this important quantity satisfies the relation 

T i = Y^w ij T s + l, (9) 

3 

where i£T and uiij is transition probability for the par- 
ticle of going from node i to node j. According to the 
definition of the random-walk problem, it is not difficult 
to know that w%j — aij/di, which is exactly the 
element of the matrix Z _1 A. 

In order to facilitate the description, we distinguish all 
nodes in M by assigning each of them a unique number. 
We label consecutively all nodes, excluding those in T, 
from 1 to N — \T\ and trap nodes are numbered from 
N - |T| + 1 to N. Then Eq. © can be rewritten in 
matrix form as 

T' = W'T' + e, (10) 

where T' = [Ti, T 2 , . . . , Tjy-|r|] T (the superscript T 
of the vector represents transpose) is an (N — 
dimensional vector, e is the (N — |r|)-dimensional unit 
vector (1, 1, . . . , 1) , and W' is the transition matrix cor- 
responding to the trapping problem. Equation (jlpp can 
be further recast as 

T' = [L']- 1 e, (11) 



In this section, we formulate the problem of random 
walks on the network M with multiple traps, which is a 
discrete-time random walk of a particle in the presence 
of several perfect traps placed on certain nodes. At each 
time step, the particle jumps with equal probability from 
its current location to one of its nearest neighbors. If the 
particle meets one of the traps, then it is absorbed. At 
last, the particle will be inevitably absorbed by the traps, 
regardless of its starting position [§0|, [Sll ■ 

It is well known that an arbitrary network can be com- 
pletely represented by its adjacency matrix. For M, its 
adjacency matrix A is a matrix consisting of entries 
or 1, with an order N x N (N is the number of nodes 
in M). The element Ojj of A is defined as fol- 

lows: fly = 1 if i and j are neighbors and a%j = 
otherwise. Then the degree di of node i is given by 
d{ = a iji t ne diagonal degree matrix Z associated 

with M is Z = diag(c?i, d 2 , ■ ■ ■ , di, . . . , d/v), and the cor- 
responding normalized Laplacian matrix of M is defined 
to be L = I — Z _1 A, where I is the identity matrix with 
order N x N. 

We use r to denote the set of traps and |T| to represent 
the number of traps. We are concerned with the expected 
time the particle spends, starting from a source node, 
before it falls on one of the traps for the first time. Let 
Ti be the expected time, frequently called first-passage 
time (FPT) or trapping time, for a particle first arriving 
at any one of the traps, given that it starts from node i. 



where 

L' = I' W' (12) 

with I' being the (N - \T\) x (N - \T\) identity matrix. 

It should be mentioned that the considered random 
walk is in fact a Markov process, and Eq. (fT2"|) is the fun- 
damental matrix of the Markov chain representing such 
an unbiased random walk. We also note that the matrix 
I' — W' on the right-hand side of Eq. (fT2"|) is actually a 
submatrix of the normalized discrete Laplacian matrix L 
of M, which is obtained from L by suppressing the last 
|r| rows and columns that correspond to the trap nodes. 

Equation (jlip shows that trapping time Tj can be ex- 
pressed in terms of the entries 1^ of inverse matrix of L' 
(i.e., a submatrix of L). Concretely, Ti is provided by 

jv— |r| 

r>= £ i^ 1 , (is) 

3=1 

which accounts for the Markov chain representing the 
random walk: The entry Z"- 1 of the fundamental matrix 
[L'] _1 for the Markov process represents the expected 
number of times the particle visits node j in the case that 
it starts off from node i (see Ref. [92] for a single trap). 
Then the MFPT (T), which is defined as the average of 
Ti over all initial nodes distributed uniformly over nodes 
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in M including the traps (93[, is given by 
iV at— |r| 

i=l i=l 

w-|r| at— |r| 

= ^ E E 

Thus, on the basis of the definition of the unbiased ran- 
dom walks, we have derived the numerical yet exact solu- 
tion to the MFPT (T) for random walks on any network 
with multiple traps, independently of the number and 
location of the traps. We note that our derivation is a 
reformulation of the backward equation satisfied by the 
MFPT and that Eq. (fl4| can also be found in the liter- 
ature in several equivalent forms [90l [9l| . 

Equation (|14p is very important, since it reduces the 
problem of calculating the MFPT (T) to computing the 
sum of the elements of the matrix [L'] _1 and can be 
used to check the results for (T) derived by other meth- 
ods, at least for networks with a small number of nodes. 
However, it is notable that although the above com- 
putational method, process, and result are applicable 
to the trapping issue on all networks, the derivation of 
Eq. (|14p requires inverting the matrix L with an order 
(N — |r|) x (N — |r|). Since the computation of inverting 
the matrix L puts heavy demands on time and mem- 
ory for large networks, by using Eq. (|14[) we can directly 
compute (T) only for networks with small size. In partic- 
ular, by applying the approach of inverting the matrix, it 
appears very difficult, even impossible, to get the exact 
dominating scaling of (T) characterizing the efficiency 
of the trapping problem. Therefore, it is of significant 
practical importance to seek an alternative method of 
computing (T) even for specific networks, which is able 
to reduce the computational effort of the method of in- 
verting the matrix. 

IV. SCALINGS OF THE MFPT FOR RANDOM 
WALKS ON MODULAR NETWORKS WITH 
MULTIPLE TRAPS 

Here we study two particular trapping problems de- 
fined in the modular scale-free networks M g . We first 
address the case that traps are located at all peripheral 
nodes; then we consider the case that traps are fixed on 
those nodes farthest from the main hub. We will show 
that the special recursive construction of the modular 
scale- free networks and the particular selections made for 
the trap locations allow for an analytical treatment of the 
MFPT to the traps. 

A. Determination of intermediate variables 

Prior to studying the MFPT to the traps, we first de- 
fine some intermediate variables and determine their val- 



ues. We denote by T g the FPT for a walker starting 
from an arbitrary peripheral node of M g to visit the hub 
for the first time and by Tf the FPT spent by a particle 
initially located at the hub to first visit any peripheral 
node. In Appendix \^ we give detailed derivations for 
Tf and Tf , which read 

„ ( 5m-2\ / m \ 9 

and 

rf = (3 ra - 8 + I^)(^) a - 2m+ 3, (16) 

respectively. 

Equations (fT5|) and (Til)]) are very useful for the follow- 
ing derivation of the exact formula for the MFPT to the 
targets. We note that Eqs. (fT5]) and ([To]) have also been 
derived in R ef. [66l| by using the technique of generat- 
ing functions [94|, but the approach used here is different 
from and relatively easier than the previous one. 

B. Exact solution to the MFPT for random walks 
with traps located at peripheral nodes 

After obtaining the intermediate quantities, we are 
now in a position to consider random walks on networks 
M g with all the (m — l) 9 peripheral nodes being occu- 
pied by traps. Our goal in this case is to determine the 
MFPT denoted by (T) g , which is the average of the FPT 
for a walker originating from a node in M g to first visit 
any target over all starting points including the traps. In 
order to find (T) g , we introduce another quantity (H) g 
defined as the FPTs of all nodes to the hub. From the 
structure of the networks, we can easily establish the fol- 
lowing recursive relations for (T} g and (H) g : 

(T) g = - ((#>,-! + T g H ) + — — — (T) g _i (17) 

TO J TO 

and 

(H) g = 1(21),-! + ^ ((TVr + If) . (18) 
Equations (TITf and (TI5)) can be rewritten as 

m(T) g - (m - l)(T) ff _! - Tf = (H) g ^ (19) 

and 

m(H) g - (#> 9 _x = (to - 1) ((T) 9 _x + Tf) . (20) 
From Eq. (|19j) . we further have 

m(T) g+1 - (to - l)(T) g - T g H +1 = (H) g , (21) 
which, together with Eqs. (Til?|) and (I2U1) . yields 

TO[TO(T) g+1 -(TO-l)(T) s -T^ 1 ] - 

[m<r> fl -(m-l)(T) ff _ 1 -T/] 
= m(H) g -{H) g _ x 

= (m-l)((T> fl _ 1+ lf) , (22) 
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that is, 



and 



(T) a+ x-(T) 



9 ™2 l"" 9 



[mT g H +1 -T g H + (m-l)T g p ). (23) 



Substituting Eqs. JTS]) and (HU) into Eq. and con- 
sidering the initial condition (T)2 = m + 1 — 2/to, we can 
solve inductively Eq. (f2"3")l to obtain the following rigorous 
expression: 



(T)„ = 



(to - l)(3m - 2) (to 2 - 2to + 2) 



TO J \ TO — 1 

2(m - l) 2 12 4 
-- i 5-^-9 -3to + 10 + ^r. 24 



Plugging this result for (T) g into Eq. (|18p and using the 
initial condition (if) 2 = to/ (to — 1), Eq. (|18p is solved to 
yield 



(H}<, 



2(3to-2)(to- l) 3 / to 



TO — 1 



2(m - l) 2 to - 1 



(5to 2 - 10m + 4)(25) 



To confirm our analytic formulas, we have compared 
them with the numerical values from the method of in- 
verting the matrix provided by Eq. (IT41 ; see Fig. [5] For 
various values of to and <?, the results for (T) g and (H) g 
obtained separately from Eqs. (|2"4")l and (|25l) are in com- 
plete agreement with those from Eq. (|14l) . This agree- 
ment serves as a mutual test of our numerical solution 
and analytical formulas, providing an important evidence 
of the validity of Eqs. (flt f . (gD and (gjj) . 

We proceed to represent (T) g and (H) g as functions of 
network size N g to obtain their dependence on N g . From 



AT 



{Ng) i-M, 



m 9 we have g = In N g / In to and m 9 /(m — l) 9 = 



-l)/lnr 



which enables us to recast Eqs. (|2"4")l 



and (|25p in terms of iV s as 



(N g )« r - 



= (iV g ) 1_ln ^ m_1 '/ ln ' 



(29) 



m _ (™- 1)(3to-2)(to 2 -2TO + 2) ^ Af u _i n(m _i)/i nm 

2(m-l) 2 lniV 12 4 

i ^.-3to + 10 + ^ (26) 

to In TO to to 



and 

{H)g = 2(3TO-2)(TO-l)3 (jVg)1 _ ln(m _ 1)/lnm 



Obviously, the exponent 77(771) is smaller than 1, show- 
ing that both (T) g and (H) g scale sublinearly with the 
network size. 

We note that the scaling in Eq. (|28p has been previ- 
ously derived in Rcf. [66] by using the theory of generat- 
ing functions, the computation process of which is a little 
complex. Equation ([28|l shows that when the hub node 
is considered an immobile trap, the trapping efficiency 
is high (even the highest among all networks [HI, l53j). 
which can be elaborated as follows. In M g the average 
distance between the hub and other nodes is only half 
of the average distance between all pairs on nodes [88| , 
suggesting that the hub node is spatially closer than any 
other node. In contrast, the degree of the hub node is 
the highest, which is why the MFPT to the hub is very 
low. 



C. Behavior of the MFPT for random walks with 
traps fixed at farthest nodes 

We now focus on the trapping problem with traps be- 
ing placed on the nodes farthest from the main hub, 
which are expected be more difficult to visit compared 
with the peripheral nodes [95| . 



1. Related definitions and quantities 



In M g the maximum value of the distance from the 
main hub to other nodes is g. We let F g denote the set of 
those nodes in M g at a distance g from the main hub of 
M g , hereafter called the farthest nodes of M g , and \F g \ 
denote the cardinality (number of elements in a set) of F g . 
By construction, M g is composed of a primal M ff _i and 

(x) 

to — 1 copies of Mg-i, denoted separately by M g _ x (x = 
1, 2, . . . , m — 1). For Mi, its farthest nodes are exactly 
its to — 1 peripheral nodes; for M2, its farthest nodes 

(x) 

correspond to the hub nodes of all M[ . Proceeding 
analogously, for g > 3 the farthest nodes of M g must 

and the farthest nodes of 



2(m-l) 2 ln7V g to-1 



m - 



In to 



(5to 2 - 10m + 4) 



(27) 



(x) 

belong to all subgraphs Mg_ v 

the primal central subgraphs (i.e., M s _2) forming M g x \ 
constitute F g . Thus, we have 

|F s | = (m-l)|F 9 _ 2 |. (30) 

Considering \F±\ — to — 1 and \Fi\ — to — 1, the recursive 
relation can be solved to obtain 



Equations ()26|) and (|27|) imply that in the limit of large 
network size (i.e., N g — > 00), both (T) g and (H) g grow 
asymptotically as power-law functions of network size N g 
with the same exponent f?(m) = 1 — ln(m — l)/lnm: 



\F n 



'(TO- 1)(9+!)/ 2 , g is Odd, 

(to — I) 5 / 2 , g is even. 



(31) 



(T) g - (Ngf^ = (7V g )l-ln(™-l)/lnm 



(28) 



Next we concentrate on the MFPT from the hub to 
the farthest nodes in M g , which will be denoted by (T)^ 
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henceforth since, as will be shown, it has the same scaling 
as that of the average of MFPTs to the farthest nodes 
F g , taken over all starting points. For a convenient de- 
scription of the computation for the MFPT to the far- 
thest nodes, we introduce more variables. For M g , let 
H g and R g express the sets of the main hub and periph- 
eral nodes, respectively. In addition, for those nodes of 

(x) 

Mg that belong to M g _ l7 we can further classify them 
in the following way. Let FL g - n (n — 1, 2, . . . , g — 1) 
denote the set of those local hubs that are directly con- 
nected to g — n classes of local peripheral nodes in P z 
and Rg- n (n — 1, 2, . . . ,g — 1) stand for the set of the 
local peripheral nodes whose neighbors are g — n differ- 
ent local hubs belonging to M z . It is easy to verify that 
the respective degrees of nodes in R g - n and H g _ n are 



K g-n = m 

respectively. 



fJ 



n and K^L n 



R g {n) and H g (n), the exact expressions for which are 



R g (n) = (m-l)"/ 2+1 



3m — 2 / m 



3m — 2 / m 



m — 1 V m — 1 



m \m — 1 
8— n— 3 

+ 1 (32) 



and 

H g (n) = (m-l)(™+ 3 )/ 2 
-(3m -2) 



3m — 2 / m 



- 2 



m — 1 



m \ m — 1 

g-n-3 

+ 2m-3, (33) 



respectively. 

Using the obtained intermediate quantities, we can de- 
rive an exact formula for (T)^. We distinguish two cases: 
(i) g is odd and (ii) g is even. For odd g we have 



2. Exact solution to the MFPT from the hub to farthest 
nodes 



According to the structure of M g , for a walker start- 
ing from the main hub, in order to reach the farthest 



nodes, it should follow the path H g 



R 



H. 



3-2 



9-3 •• • 



-> R 



g-(n-2) 



-> H, 



Rg 

S-(n-l) 



Hg-1 
Ra—n 



g-(n+l) 



R g _r n+ 2) ' ' • Ri(otHi). Then it is natural 
to define the following quantities. Let R g {n) and H g {n) 
represent, respectively, the FPT from a node in R g - n to 
any of its neighboring nodes in i? g _(„ + i) and the FPT 
from a node in Ff g - n to any of its neighbors belonging to 



( ff -l)/2-l (g-l)/2-2 
i=0 i=0 



+ (m-l)i? g (.g-3) + m 



(34) 



By plugging Eqs. ([32]) and (|33J) into Eq. dMJ) and do- 
ing some algebra, we find a closed-form solution to (T)^ 
given by 



(T)" = 



3m - 2 



( m \ J i -,N(g+i)/2 2m / 
7 (m-l) WT " -(m-1) 

(m - l)(3m - 2) (m 2 - 2r?x + 3) / m 



(s+i)/2 



(m- l)g + 



(m - 2) (2m - 1) \m — 1 

(m - 1) (3m 3 - 9m 2 + 14m - 4) 



(m - 2) (2m - 1) 



(35) 



R 



g-(n+l)- 



In Appendix [B] we report the derivation for When g is even, it is not difficult to reach the following 
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FIG. 3: (Color online) Mean first-passage time (T)f as a 
function of generation jona log-log scale for two special cases 
of m = 3 and 4. The open symbols indicate the numerical 
results obtained by direct calculation from Eq. (|13|) ; the solid 
symbols display the analytical values provided by Eqs. (|35p 
and ((37)) . 



expression: 

9/2-2 g/2-2 

< T >f = T a H + E R 9 w+ E ^(2*+!) 

i=0 i=0 

+H g {g-3)+m. (36) 

Inserting Eqs. (|32p and (|33j) into Eq. after some al- 
gebra, the explicit expression for (T)^ is obtained, which 
reads 



(T)" = 



2(3m - 2) / m 



(m - l) s 



(m — 2)m \m - 1 
(m - l)(3m - 2) (m 2 - 2m + 3) 



-(m-1) 



9/2+1 



(m - l)g + 



(m - 2)(2m - 1) 

3m 4 - 11m 3 + 19m 2 



14m + 4 



(m - 2)(2m - 1) 



(37) 



To check the validity of Eqs. ([33]) and ([3T|). we also com- 
pute (T 1 )^ numerically by using the approach of inverting 
the related matrix; see Eq. (TIB")) . The results obtained by 
analytical and numerical methods completely agree with 
each other. The comparison is shown in Fig. [3] Equa- 
tion (1351) . together with Eq. (|37|) . indicates that for large 
networks, i.e., N g —> oo, 

(T)f - (^ 3 ) S(m) = (iV 9 )l-ln(™-l)/(21nm) ; (33) 

with the exponent 9(m) = 1 — ln(rn — l)/(2 lnm) smaller 
than 1. 

Thus far we have found the rigorous formula for the 
MFPT (T}£ to farthest nodes in M g and its dependence 
on network size N g . We stress that the analytical com- 
putation for the MFPT (T) g to the farthest nodes that 
average all starting points in M g is rather lengthy and 



awkward. However, it is easy to infer that when g is 
large enough, the dominant term of (T) g also increases 
as a power-law function of network size N g with an expo- 
nent identical to that of (T)^ , which can be understood 
from the following heuristic explanation. Note that M g 
consists of to subgraphs, which are copies of M s _i. For 
those nodes in the central subgraph, their MFPT to the 
farthest nodes is equal to (T)^ + (if) 9 _i, the dominant 
term of which is (T)^; for nodes in each of the m — 1 
fringe subgraphs (a: = 1, 2, . . . , m-1), their MFPT 

to the farthest nodes is identical but smaller than (T)^. 
Hence, for all nodes in M g , the dominating term of the 
MFPT (T) g is proportional to (A^M but its prefactor 
may be different from that of (T)^. 



D. Result analysis 



Equations (f2"5|) and ([35]) show that when traps are posi- 
tioned at several particular nodes, the MFPTs to the tar- 
get node are very small, which scale sublinearly with the 
network order. When either peripheral nodes or farthest 
nodes are occupied by traps, the characteristic exponent 
rj(m) or 9{m) is a decreasing function of m: When the 
parameter to increases from 3 to oo, both r/(m) and 9(m) 
drop and are close to zero. Therefore, the efficiency of the 
random-walk process is reliant on m: The larger the pa- 
rameter m, the more efficient the random-walk process. 
The fact that both trapping processes are very efficient 
demonstrates that the modular scale-free networks be- 
ing studied exhibit an efficient configuration for random 
walks with traps positioned at certain given nodes. 

In contrast, for each given parameter to, rj(m) is 
smaller than 9(m), which implies that when traps are 
located at peripheral nodes, the trapping efficiency is 
higher than that of the case when traps are placed on 
farthest nodes. Thus, the two trapping processes defined 
on the networks under consideration display rich behav- 
ior in the context of MFPTs to the traps. The difference 
between ?y(m) and 6(m) shows that the number and lo- 
cation of traps sensitively affect the behavior of random 
walks on the modular scale-free networks. 

Actually, the intrinsic structure of the modular scale- 
free networks is responsible for the high efficiency of 
random walks performing on them with certain nodes 
being occupied by traps. In these networks, there are 
many small highly integrated clusters, which group into 
a few larger but less compact modules linked by local 
hub nodes; see Fig. [TJ These relatively large modules 
combine to form even larger and fewer groups, which are 
further joined to shape a fine modular and scale-free ar- 
chitecture, a topology that accounts for the fast diffusion 
phenomenon in M g . 

In the case that traps are placed on peripheral nodes, 
when a particle originates from a node in a duplicate 
M g -i (an element of M g ), it will either be directly 
trapped by one of the traps or jump to local hub nodes in 
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a few steps. These local hubs play a bridge role linking 
different small modules together. After arriving at local 
hub nodes, the particle can be easily trapped in a short 
time. In contrast, if the particle starts off from a node 
in the original M 9 _i (the central part of M g ), it will eas- 
ily visit local hub nodes or the hub first, through which 
it can find the way to one of the traps quickly. Thus, 
the particle can drop into the traps very fast wherever it 
starts to jump, which can be understood from the above 
heuristic argument based on the inherent structure of the 
considered networks. 

When the traps are fixed on the farthest nodes, to find 
a garget, the walker must first visit the local hubs and lo- 
cal peripheral nodes of a larger and sparser cluster, start- 
ing from which it continues to arrive at the local hubs 
and local peripheral nodes of smaller and denser groups. 
From Eqs. (|3"2"j) and (|3U)l we know that the expected time 
between local hub nodes and local peripheral nodes in in- 
ner subgraphs rely on their size or deepness (i.e., g — n): 
The smaller the value of g — n, the smaller the size of 
inner subgraphs, and the higher the expected time. This 
can account for the main reason the farthest nodes are 
more difficult to reach than the peripheral nodes. 



very specific sets of traps and is hard to generalize to 
other sets of traps. 
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Appendix A: Derivation of Tg and 

According to the particular structure of the networks, 
for any g > 1, the two quantities Tg and Tjf obey the 
following recursion relations: 



T p = 

9 



1 



(m - 2) + g 
l + T H 



9-1 

E 



1 + (to - 2) 

^9 P ) 



1 + Ti 



(Al) 



and 



V. CONCLUSIONS 

We have studied the random- walk dynamics on a fam- 
ily of modular scale-free networks with multiple traps, 
which exhibit remarkable characteristics observed for var- 
ious real-life networks, such as social and biological net- 
works. We first deduced a general formula for the MFPT 
to the traps in a generic network, which is expressed in 
terms of several elements of a matrix associated with the 
trapping problem. Then we studied the MFPT for two 
trapping issues on the studied networks with two differ- 
ent arrangements of targets. In the first case, peripheral 
nodes are treated as traps; in the second case, farthest 
nodes work as traps. 

For the two trapping problems, we studied both numer- 
ically and analytically the MFPT to traps, the results 
of which are compatible with each other. Our results 
show that in both cases, the MFPT varies as a power- 
law function of network size with the exponent depend- 
ing on the parameter to, which is lower than 1 in the 
full range of to. Thus, the studied networks display an 
efficient architecture in favor of diffusion. Moreover, we 
demonstrated that, compared with the second case, the 
diffusion is faster in the first case, which indicates that 
the transport efficiency relies on the number and loca- 
tion of the absorbing nodes. We also showed that the 
modular topology, together with the scale-free behavior, 
is responsible for the quick diffusion processes, as well as 
the scaling difference of the MFPT for the two trappings 
running on the networks addressed. We expect that our 
work can provide insight into designing networks with a 
structure in favor of diffusion. Finally, it should be men- 
tioned that the method developed here applies only to 



T, 



H 



Eli (m-l) 

9-1 



(m-l) 9 + 



i=l 



iy (1 + if + T t 



(A2) 



The three terms on the right-hand side (rhs) of 
Eq. (jAip can be explained as follows. The first term 
is based on the fact that the walker takes only one time 
step to first reach the hub. The second term describes 
the process by which the particle first jumps to one of its 
to — 2 neighbors belonging to P in one time step and then 
takes more steps to first get to the target node. The 
last term accounts for the fact that the walker first makes 
a jump to a local hub node belonging to H 2 , then takes 
T? 1 time steps, starting off from the local hub, to reach 
any node in P, and continues to jump Tg more steps to 
reach the target node for the first time. 

Analogously, the two terms on the rhs of Eq. (IA2[) are 
based on the following two processes. The first term de- 
scribes the fact that the walker, starting from the hub, 
requires only one time step to hit a peripheral node. 
The second term explains such a process that the walker, 
starting off from the hub, first jumps to a local peripheral 
node belonging to ¥ z in one time step, then makes Tf 
jumps to the hub, and proceeds to any node in P, taking 
T^ more time steps. 

After merging similar items, Eqs. (|A1[) and (IA2[) can 
be rewritten as 



T P = (m-2) 



9-1 

1=1 



(A3) 



9 



and 



1 



3-1 



^(m-lf + ^im-lfTf 



i=l 



j,H _ 

9 (m - l)s 



respectively. Equations (|A3|) and (|A4I) lead to 
T Q P +1 -T Q P = 1 + T Q H 



and 



rpH 

1 9+l 



1 r g = i + 1 r p . 

m — 1 9 m — 1 9 



(A4) 
(A5) 

(A6) 



According to Eq. (|A6I) . we obtain 
1 



jt m — 1 a 



+9+1 +9+1 



x Jig 

TO — 1 3 



(A7) 



where the relation provided in Eq. (|A5|) was used. Apply- 

-2 



ing the initial condition = ^rr, we solve Eq. (|AT|) 



to obtain 



5m - 2 



TO — 1 



- 1. 



(A8) 



Inserting the result for into Eq. (|A5|i . we arrive the 



exact formula for T p given by 
7m - 2 



Ti 



3m -8 + 



m — 1 



2m + 3. (A9) 



Appendix B: Derivation of R g (n) and H g (n) 



of the m — 2 neighboring nodes belonging to R g - n , from 
which it continues to jump R g (n) steps to first visit a 
target; this is accounted for by the first term on the rhs. 
Alternatively, the walker can go to a local hub belong- 
ing to H g _/ n _i\, then takes time H g (n — 1) to reach a 
neighbor in R g - n , and proceeds to bounce R g {n) steps 
to hit a target for the first time, this process is explained 
by the second term. The third term describes the pro- 
cess by which the walker goes directly to a target node. 
The fourth term represents the process that the walker 
first jumps to a neighbor belonging to H g _/ g _i), makes a 
move returning to a node in R g - n , and then walks con- 
tinuously in time R g (n) to arrive at a destination node. 
Finally, the last sum term explains the fact that the par- 
ticle goes to a local hub in H g _i (2 < i < g — (n + 2)) 
from which it takes an average time T* 1 to return to 
one of its neighbors in R g - n , and then moves on aver- 
age R g -n steps to get to a target. Analogously, we can 
explain Eq. (|B"2|) . 

After some algebra, Eqs. (|F31[) and ()B2|) can be simpli- 
fied to 

g-{n+2) 

R g (n) =m-2+g-n+H g {n-l) + l+ ^ T? 1 , (B3) 

i=2 



and 



H g (n) = (m-l)' i+1 - 9 ^(m-l) i + (m-l)i? s (n-l) 

i=l 

+ (m - l) n + 3 s + 

g-(n+2) 

( m - Yl [(m-iyTf}. (B4) 



For the two quantities R g (n) and H g (n), the following 
relations hold: 



R g {n) 



and 



1 



K R 

g—n 



{(m-2)[l + R g (n)} + 



[1 + H g (n - 1) + R g (n)} + 1 + [2 + R g {n)] 

S-O+2) 

+ ll + T? + R g (n)}} (Bl) 



Inserting Eq. (IB4[) into Eq. (|B3I) and utilizing the ini- 
tial condition R g (Q) = m - 1 + g + T g ff + Y%L% T" = 

(m-l)(3m-2)(m 2 -m+l) / rn 
m \ m — 1 

solved to get 



2to + 3, Eq. (|B4l is 



R g (n) = (m-l)"/ 2+1 



3m — 2 / m 



3m — 2 / m 



to — 1 V m — 1 



m \m — 1 



- 2 



(B5) 



^» = ] ^{(m-l) 9 -"[l + J R 3 (n-l) + // s (n)] 



(to - + (m - l)[m + Hg(n)] 

g-(n+2) 

+ Y (m-mi + Tf + Hg(n)}}. (B2) 



i=2 



Equation (|B1[) can be elaborated as follows. Originat- 
ing from a node in R g - n , the particle can jump to one 



Substituting this expression for R g (n) into Eq. (|B4j) and 
solving Eq. (|B4|) . we obtain 



H g (n) = (m-1) 



(n+3)/2 



3m — If m 



-(3m- 2) 



rn — 1 



m \m — 1^ 
9— n-3 

+ 2m - 3 . (B6) 
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